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Conformal mapping has been applied mostly to harmonic functions, i.e. solu- 
tions of Laplace's equation. In this paper, it is noted that some other equations 
are also conformally invariant and thus equally well suited for conformal map- 
ping in two dimensions. In physics, these include steady states of various nonlin- 
ear diffusion equations, the advection-diffusion equations for potential flows, and 
the Nernst-Planck equations for bulk electrochemical transport. Exact solutions 
for complicated geometries are obtained by conformal mapping to simple geome- 
tries in the usual way. Novel examples include nonlinear advection-diffusion layers 
around absorbing objects and concentration polarizations in electrochemical cells. 
Although some of these results could be obtained by other methods, such as Boussi- 
nesq's streamline coordinates, the present approach is based on a simple unifying 
principle of more general applicability. It reveals a basic geometrical equivalence of 
similarity solutions for a broad class of transport processes and paves the way for 
new applications of conformal mapping, e.g. to non-Laplacian fractal growth. 
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1. Introduction 

Complex analysis is one of the most beautiful subjects in mathematics, and, in 
spite of involving imaginary numbers, it has remarkable relevance for 'real' appli- 
cations. One of its most useful techniques is conformal mapping, which transforms 
planar domains according to analytic functions, w = f{z), with f'{z) ^ 0. Geomet- 
rically, such mappings induce upon the plane a uniform, local stretching by |/'(2;)| 
and a rotation by arg f'{z). This 'ampli-twist' interpretation of the derivative im- 
plies conformality, the preservation of angles between intersecting curves (Needham 
1997). 

The classical application of conformal mapping is to solve Laplace's equation, 

V^^-O, (1.1) 
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i.e. to determine harmonic functions, in complicated planar domains by mapping 
to simple domains. The method relies on the conformal invariance of Eq. (1.1), 
which remains the same after a conformal change of variables. Before the advent of 
computers, important analytical solutions were thus obtained for electric fields in 
capacitors, thermal fluxes around pipes, inviscid flows past airfoils, etc. (Needham 
1997; Churchill & Brown 1990; Batchelor 1967). Today, conformal mapping is still 
used extensively in numerical methods (Trefethen 1986). 

Currently in physics, a veritable renaissance in conformal mapping is centering 
aroimd 'Laplacian-growth' phenomena, in which the motion of a free boimdary is 
determined by the normal derivative of a harmonic function. Continuous problems 
of this type include viscous fingering, where the pressure is harmonic (Saffman & 
Taylor 1958; Bensimon et al. 1986; Saffman 1986), and solidification from a super- 
cooled melt, where the temperature is harmonic in some approximations (Kessler 
et al. 1988; Cummings et al. 1999). Such problems can be elegantly formulated 
in terms of time-dependent conformal maps, which generate the moving boundary 
from its initial position. This idea was first developed by Polubarinova-Kochina 
(1945a; 1945b) and Galin (1945) with recent interest stimulated by Shraiman & 
Bensimon (1984) focusing on finite-time singularities and pattern selection (How- 
ison 1986; Tanveer 1987; Dai et al. 1991; Ben Amar 1991; Howison 1992; Tanveer 
1993; Ben Amar & Brener 1996; Ben Amar & Poire 1999; Feigenbaum et al. 2001). 

Stochastic problems of a similar type include diffusion-limited aggregation (DLA) 
(Witten & Sander 1981) and dielectric breakdown (Niemeyer et al. 1984). Recently, 
Hastings & Levitov (1998) proposed an analogous method to describe DLA using 
iterated conformal maps, which initiated a flurry of activity applying conformal 
mapping to Laplacian fractal-growth phenomena (Davidovitch et al. 1999, 2000; 
Barra et al. 2002a, 2002b; Stepanov & Levitov 2001; Hastings 2001; Somfai et al. 
1999; Ball & Somfai 2002). One of our motivations here is to extend such pow- 
erful analytical methods to fractal growth phenomena limited by non-Laplacian 
transport processes. 

Compared to the vast literature on conformal mapping for Laplace's equation, 
the technique has scarcely been applied to any other equations. The difficulty with 
non-harmonic functions is illustrated by Helmholtz's equation, 

V20 = 0, (1.2) 

which arises in transient diffusion and electromagnetic radiation (Morse & Feshbach 
1953). After conformal mapping, w = f{z), it acquires a cumbersome, non-constant 
coefficient (the Jacobian of the map), 

l/?VV = .^. (1.3) 

Similarly, the bi-harmonic equation, 

V^V^ = 0, (1.4) 

which arises in two-dimensional viscous flows (Batchelor 1967) and elasticity (Muskhe- 
lishvili 1953), transforms with an extra Laplacian term (see below), 

|/'|4v2 V^,/. = -4 If'fV^cp. (1.5) 
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In this special case, conformal mapping is commonly used (e.g. Chan et al. 1997; 
Crowdy 1999, 2002; Barra et al. 2002b) because solutions can be expressed in terms 
of analytic functions in Goursat form (Muskhclishvili 1953). Nevertheless, given the 
singular ease of applying conformal mapping to Laplace's equation, it is natural to 
ask whether any other equations share its conformal invariance, which is widely 
believed to be unique. 

In this paper, we show that certain systems of nonlinear equations, with non- 
harmonic solutions, are also conformally invariant. In section 2, we give a simple 
proof of this fact and some of its consequences. In section 3, wc discuss applications 
to nonlinear diffusion phenomena and show that single conformally invariant equa- 
tions can always be reduced to Laplace's equation (which is not true for coupled 
systems). In section 4, we apply conformal mapping to nonlinear advection-diffusion 
in a potential flow, which is equivalent to streamline coordinates in a special case 
(Boussinesq 1905). In section 5, we apply conformal mapping to nonlinear electro- 
chemical transport, apparently for the first time. In section 6, we summarize the 
main results. Applications to non-Laplacian fractal growth are a common theme 
throughout the paper. (See sections 2, 4, and 6.) 



2. Mathematical Theory 

(a) Conformal Mapping without Laplace's Equation? 
The standard application of conformal mapping is based on two facts: 

1. Any harmonic function, (j). in a singly connected planar domain, is the 
real part of a analytic function, the 'complex potential' (which is unique 
up to an additive constant): = Re $(w). 

2. Since analyticity is preserved under composition, harmonicity is preserved 
under conformal mapping, so = Re '^{f{z)) is harmonic in 0.^ = /~^(0^). 

Presented like this, it seems that conformal mapping is closely tied to harmonic 
functions, but Fact 2 simply expresses the conformal invariance of Laplace's equa- 
tion: A solution, (f){w), is the same in any mapped coordinate system, <f){f{z)). Fact 
1, a special relation between harmonic functions and analytic fimctions, is not really 
needed. If another equation were also conformally invariant, then its non-harmonic 
solutions, (j){w,w), would be preserved under conformal mapping in the same way, 
</,(/(z),7(I)). (Sec Fig. 1.) 

In order to seek such non-Laplacian invariant equations, we review the transfor- 
mation properties of some basic differential operators. Following Argand and Gauss, 
it is convenient to represent two-dimensional vectors, a = a^x + ayy, as complex 
numbers, a = Qx + ayi. We thus express the gradient vector operator in the plane 
as a complex scalar operatorf , 

^ ^ d ^ d „ d d 

V = x—+y— ^ v = — +z — , 2.1 
ox ay ox oy 

which has the essential property that V/ = if and only if / is analytic, in which 
case, V/ = 2/' (Needham 1997). Since a • b = Re afo, the Laplacian operator can 

t Although ^ = is more common in the mathematical literature, we prefer V for appli- 
cations in transport theory because gradients play a central role. 
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<^(w,w) (\>(f(z)J(z)) 

Figure 1. Conformal mapping, w = f{z), of a solution, 4>, to a conformally invariant 
equation from a complicated domain, fi^, and a simple domain, fi^. 



be expressed as, V • V = Rc VV = VV (if mixed partial derivatives can be taken 

in any order). Similarly, the 'advection operator', which acts on two real functions 

^ and c, takes the form, V<^ • Vc = Re (V(^)Vc. 

Under a conformal mapping of the plane, w = f{z), the gradient transforms as, 
= /' Vu,. This basic fact, combining the ampli-twist property and the chain rule, 

makes it easy to transform differential operators (Needham 1997). The Laplacian 

transforms as, 

V,V, = (V,/')V^ + |/?V^V^ = l/f V^V,„ (2.2) 

where Vz/' = because /' is also analytic. This immediately implies the confor- 
mal invariance of Laplace's equation (1.1), and the non-invariance of Helmholtz's 
equation (1.2). The transformation of the bi-harmonic equation (1.4) in Eq. (1.5) 
is also easily derived with the help of Needham's identity, A|/p = 4|/'p, applied 
to /'. 

Everything in this paper follows from the simple observation that the advection 
operator transforms just like the Laplacian, 

Re (Vz(/.) Vzc = l/f Re (V^0) V^c. (2.3) 

Each operator involves a 'dot product of two gradients', so the same Jacobian fac- 
tor, l/'P, appears in both cases. The transformation laws, Eq. (2.2) and Eq. (2.3), 
are surely well known, but it seems that some general implications have been over- 
looked, or at least not fully exploited in physical applications. 



(6) Conformally Invariant Systems of Equations 

The identities (2.2) and (2.3) imply the conformal invariance of any system of 
equations of the general form, 

N / N \ 

^ ( ai{(l>) V^ci>i + ^ ai,(</)) V4>i ■ V4>j 1=0 (2.4) 

where the coefficients aj(0) and aij{4>) may be nonlinear functions of the unknowns, 
<p = {^i,(f)2, . . . ,(Pn), but not of the independent variables or any derivatives of the 
unknowns. Thus we arrive at our main result: 
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Theorem 2.1. (Conformal Mapping Theorem.) Let (p{w,W) satisfy Eq. (2.4) in 
a domain Slw of the complex plane, and let w = f{z) be a conformal mapping from 
flz to flyj. Then (j)[f{z),f{z)) satisfies Eq. (2.4) in fl^. 

Whenever the system (2.4) can be solved analytically in some simple domain, 
the Theorem produces a family of exact solutions for all topologically equivalent 
domains. Otherwise, it allows a convenient numerical solution to be mapped to 
more complicated domains of interest. This is an enormous simplification for free 
boundary problems, where the solution in an evolving domain can be obtained by 
time-dependent conformal mapping to a single, static domain. 

Conformal mapping is most useful when the boundary conditions are also invari- 
ant. Dirichlet {(j^i = constant) or Neumann (n • V^i = 0) conditions are typically 
assumed, but here we consider the straight-forward generalizations, 

JV 

6i(0) =0 and ^ 6^(0) (n • = (2.5) 

respectively, where bi{<p) and bij{<p) are nonlinear functions of the unknowns, aj 
is a constant, and ft is the unit normal. The conformal invariance of the former is 
obvious, so we briefly consider the latter. 

It is convenient to locally transform a vector field, F, along a given contour 
as, F = tF, so that Re F and Im F are the projections onto the unit tangent, 
t = dz/\dz\, and the (right-handed) imit normal, n = —it, respectively. Since the 
tangent transforms as, tw = tzf'/\f'\, and the gradient as, Vz = /'Vu.. wc find, 

= \f'\Vw The invariance of Eq. (2.5) follows after taking the imaginary part 
on the boundary contour. 

(c) Gradient-Driven Flux Densities 

Generalizing V(/) for Laplacian problems, we define a 'flux density' for solutions 
of Eq. (2.4) to be any quasi-linear combination of gradients, 

N 

Fi = ^Ci,(0)V<^,-, (2.6) 

i=i 

where Cij{(p) are nonlinear functions. The transformation rules above for the gra- 
dient apply more generally to any flux density, 

F,=JF^ and F, = \f'\F^. (2.7) 

These basic identities imply a curious geometrical equivalence between solutions to 
different conformally invariant systems: 

Theorem 2.2. (Equivalence Theorem.) Let <p^^^ and cf)'^^^ satisfy equations of the 
form (2.4) with corresponding flux densities. F^^^ and F^'^\ of the form (2.6). If 
Fz^^ = aF^^^ on a contour Cz for some complex constant a, then F^^ = aF^^ on 
the image, Cyj = f{Cz), after a conformal mapping, w = f{z). 

An important corollary pertains to 'similarity solutions' of Eqs. (2.4) and (2.5) in 
which certain variables {(pi} involved in a flux density depend on only one Cartesian 
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coordinate, say Hew, after conformal mapping: (pi = Gi(Re f{z)). (Our examples 
below are mostly of this type.) Such special solutions share the same flux lines 
(level curves of Im f{z)) and iso-potcntials (level curves of Rc f{z)) in any geometry 
attainable by conformal mapping. They also share the same spatial distribution of 
flux density on an iso-potential, although the magnitudes generally differ. 

An important physical quantity is the total normal flux through a contour, often 
called the 'Nusselt number', Nu . For any contour, C, we define a complex total flux, 
7(C) = J^F \dz\ = J^Fdz, such that ReI{C) is the integrated tangential flux and 
Im/(C) = Nu(C). From Eq. (2.7) and dw = f'dz, wc conclude, I{C,) = I{Cu,). 
Therefore, flux integrals can be calculated in any convenient geometry. 

This basic fact has many applications. For example, if is constant on a 
contour Cw ~ f{Cz), then for any conformal mapping, wc have, I{Cz) = I{Cw) = 
i{Cyj)Fyj, where i{Cyj) is simply the length of Cyj. For fluxes driven by gradients 
of harmonic functions, this is the basis for the method of iterated conformal maps 
for DLA, in which the 'harmonic measure' for random growth events on a fractal 
cluster is replaced by a uniform probability measure on the unit circle (Hastings & 
Levitov 1998). 

More generally, a non-harmonic probability measure for fractal growth can be 
constructed for any flux law of the form (2.6) for flelds satisfying Eq. (2.4). Ac- 
cording to the results above, the growth probability is simply proportional to the 
normal flux density on the unit circle for the same transport problem after confor- 
mal mapping to the exterior of the unit disk. A nontrivial example is given below 
in section 4(c). This allows the Hastings-Levitov method to be extended to a broad 
class of non-Laplacian fractal-growth processes (Bazant, Choi & Davidovitch 2003). 



The Conformal Mapping Theorem is even more general than it might appear 
from our proof: The domain, Clz, may be contained in any two-dimensional manifold. 
This becomes clear from the recent work of Entov and Etingof (1991; 1997), who 
solved viscous fingering problems on various curved surfaces by conformal mapping 
to the complex plane, e.g. via stereographic projection from the Riemann sphere. 
They exploited the fact that Laplace's equation is invariant under any conformal 
mapping, w = f{z), from the plane to a curved surface because the Laplacian 
transforms as = J V^, where J{f{z)) is the Jacobian. The system (2.4) shares 
this general conformal invariance because the advection operator transforms in the 
same way, V^(/)- V^c = J V^0- V^c. The application of these ideas to non-Laplacian 
transport-limited growth phenomena on curved surfaces is work in progress with J. 
Choi and D. Crowdy; here we focus on conformal mappings in the plane, described 
by analytic functions. 

3. Physical Applications to Diffusion Phenomena 

Conformally invariant boundary-value problems of the form (2.4) and (2.5) com- 
monly arise in physics from steady conservation laws, 



(d) Conformal Mapping to Curved Surfaces 








(3.1) 
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for gradient-driven flux densities, Eq. (2.6), with algebraic {b{ci) = 0) or zero-flux 
(n • Fj = 0) boundary conditions, where c, is the concentration and F, the flux of 
substance i. Hereafter, we focus on flux densities of the form, 

Fj = Cj Uj - Di{ci) Vci, Uj cx (3.2) 

where Di{ci) is a nonlinear diffusivity, Uj is an irrotational vector field causing ad- 
vcction, and (p is a (possibly non-harmonic) potential. Examples include s advection- 
diffusion in potential flows and bulk electrochemical transport. 

Before discussing these cases of coupled dependent variables, it is instructive to 
consider nonlinear difFusion in only one variable. The most general equation of the 
type (2.4) for one variable is, 

a(c) V^c= |Vc|2. (3.3) 

This equation arises in the Stefan problem of dendritic solidification, where c is the 
dimensionless temperature of a supercooled melt and a(c) is Ivantsov's function, 
which implicitly determines the position of the liquid-solid interface via a(c) = 1 
(Ivantsov 1947). In two dimensions, Bedia & Ben Amar (1994) prove the confor- 
mal invariance of Eq. (3.3) and then study similarity solutions, c(^,?7) = G(r?), by 
conformal mapping, w = ^ + i-q, to a, plane of parallel flux lines, 

a{G)G" = {G'f, (3.4) 

where an ordinary differential equation is solved. 

More generally, reversing these steps, it is straight-forward to show that any 
monotonic solution of Eq. (3.4) produces a nonlinear transformation, c = G{(p), 
from Eq. (3.3) to Laplace's equation (1.1), which implies conformal invariance. 
There are several famous examples. For steady concentration-dependent diffusion, 

V • {D{c)Vc) = 0, (3.5) 

it is Kirchhoff's transformation (Crank 1975), (p = G~^{c) = D{x)dx. For Burg- 
ers' equation in an irrotational flow (u = — V/i), 

-I- Au • Vu = uV'^u, (3.6) 

at 

which is equivalent to the KPZ equation without noise (Kardar et al. 1986), 

^=uV^h+^\Vhf, (3.7) 

it is the Cole-Hopf transformation (Whitham 1974), (j) = G-'^{h) = e^^/'^" , which 
yields the diffusion equation, ^ = u'V'^cj), and thus Laplace's equation in steady 
state. 

In summary, the general solutions to Equation (3.3) are simply nonlinear func- 
tions of harmonic fimctions, so, in the case of one variable, our theorems can be 
easily understood in terms of standard conformal mapping. For two or more coupled 
variables, however, this is no longer true, except for special similarity solutions. The 
following sections discuss some truly non-Laplacian physical problems. 



Article to appear in Proc. Roy. Soc. A (submitted on 15 January 2003). 



8 



M. Z. Bazant 



4. Steady Advection-Diffusion in a Potential Flow 

We begin with a well known system of the form (2.5), the only one to which con- 
formal mapping has previously been applied (see below), albeit not in the present, 
more general context. Consider the steady diffusion of particles or heat passively 
advected in a potential flow, allowing for a concentration-dependent diffusivity. For 
a characteristic length, L, speed, U, concentration, C, and diffusivity, D{C), the 
dimensionless equations are 

VV = and Pe V?!> • Vc = V • (6(c) Vc), (4.1) 

where (p is the velocity potential (scaled to UL), c is the concentration (scaled to 
C), 6(c) is the dimensionless diffusivity, and Pe — UL/D is the Pcclet number. The 
latter equation is a steady conservation law for the dimensionless flux density, F = 
Pe cV(p — b{c)'Vc (scaled to DC/L). For 6(c) = 1, these classical equations have 
been studied recently in two dimensions, e.g. in the contexts of tracer dispersion in 
porous media (Koplik et al. 1994, 1995), vorticity diffusion in strained wakes (Fames 
& Bush 1999; Hunt & Fames 2002), thermal advection-diffusion (Morega & Behan 
1994; Sen & Yang 2000), and dendritic solidification in flowing melts (Kornev & 
MukhamaduUina 1994; Cummings et al. 1999). 

(a) Similarity Solutions for Absorbing Leading Edges 

Let us rederive a classical solution in the upper half plane, w = £, + irj (77 > 0) , 
which we will then map to other geometries. As shown in the top left panel of Fig. 2, 
consider a straining velocity field, ^ = Re$, $ = w^, u = $' = 2W = 2^ — 2ir], 
which advects a concentrated fluid, c(^, 00) = 1, toward an absorbing wall on the 
real axis, c(^, 0) = 0. Since the ry-component of the velocity (toward the wall) is 
independent of ^, as are the boundary conditions, the concentration depends only 
on T]. The scaling function, c(^, ry;Pe) = S{-\/Per]) = S{fi), satisfies 

-2fiS' = {b{S) Sy, S{0) = 0, 5(00) = 1, (4.2) 

which is straightforward to solve, at least numerically. For 6(5') = 1, Equation (4.2) 
has a simple, analytical solution, S{fi) = erf {fj) (e.g. Cummings et al. 1999). 

If extended to the entire ui-plane, where two fluids of different concentrations 
flow towards each other, this solution also describes a Burgcas' vortex sheet under 
uniform strain (Burgers 1948). In that case, ((/)^ c) is a three-dimensional ve- 
locity held satisfying the Navier-Stokes equations, and Pe is the Reynolds number. 
Inserting a boundary, such as the stationary wall on the real axis, however, is not 
consistent with Burgers' solution because the no-slip condition cannot be satisfied. 
The wall is crucial for conformal mapping to other geometries because it enables 
singularities to be placed in the lower half plane. 

For every conformal map to the upper half plane, w = f{z), we obtain a solution, 

(j) = Refizf and c=S (VPelm f{z)^ for Im f{z) > (4.3) 

which describes the nonlinear advection-diffusion layer in a potential flow of concen- 
trated fluid around the leading edge of an absorbing object. For a linear diffusivity, 
S{fi) = erffi, various examples are shown in Fig. 2. The choice, f{z) = \fz — a, in 
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Figure 2. Concentration profiles (contour plots) and potential-flow streamlines (yellow) 
for steady, linear advection-diffusion layers around various absorbing surfaces (gray) at 
Pe — 1. All solutions are given by Eq. (4.3), where w = f{z) is a conformal map to the 
upper half plane (top left). The color scale applies to all panels in Figs. 2 and 3. 
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Figure 3. The steady linear advection-diffusion layer around a cylindrical rim on a flat 
plate at Pe = 0.1 (left) and Pe = 10 (right). 



the middle left panel, describes a parabolic leading edge, x = {y/2a)'^ — , where 
a = Ima > 0. The limit of uniform flow past a half plate (a = 0), in the upper 
right panel is a special case discussed below. 

Another classical map, f{z) = z^/'^'^'"~f^\ describes a wedge of opening angle, 
/3, as shown in the middle right panel for (3 = tt/2 (after a rotation by 7r/4). The 
half plate {(3 = 0) and the flat wall (/3 ~ tt) discussed above are special cases. The 
diffusive flux on the surface from Eq. (4.8), |V(/)| (x \/Pe"r~'', is singular for acute 
angles, /3 < tt. The geometry-dependent exponent, v = (tt — /3)/(27r — /3), is the same 
for pure diffusion to the wedge, (j)d oc Im/(z) (Barenblatt 1995). This insensitivity 
to Pe is a signature of the Equivalence Theorem, as explained below. 

The less familiar mapping, f{z) = z^l"^ + which plays a crucial role in 

non-Laplacian growth problems (see below) , places a cylindrical rim on the end of a 
semi-infinite flat plate, as shown in the lower left panel. The solution has a pleasing 
form in polar coordinates. 



1 



r H — I cos Q 

r . 



erf 



(4.4) 
(4.5) 



where we have shifted the velocity potential, <i> = /(z)^ — 2 = z 



^. Far from 



the rim, we recover the half-plate similarity solution, since f{z) ~ ^/z as \z\ oo, 
but close to the rim, as shown in Fig. 3, there is a nontrivial dependence on Pe . 
For Pe ^ 1, a boundary layer of 0(Pc ^^/^) thickness forms on the front of the rim 
and extends to within an 0(Pc ^^/^) distance from the rear stagnation point. 

The flux density is easily calculated in the w-plane and then mapped to the 
z-plane using Eq. (2.7): 



= 2/'(z)/(z)Pe5 \/Pelm/(z) - /'(z) VJ^S*' VPe Im /(z 



(4.6) 
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where the first term describes advection and the second, diffusion. The Unes of 
advective and diffusive flux, which are level curves of Ini/(^!)^ and Re/(2:), respec- 
tively, arc independent of Pc and 6(c), as required by the Equivalence Theorem. 
In particular, the diffusive flux lines have the same shape for any flow speed or 
nonlinear diffusivity as in the case of simple linear diffusion (Pe = 0, 6(c) = 1, 
c oc Im/(z)), even though advection and nonlincarity affect the lines of total flux. 

The lines of total flux, called 'heatlines' in thermal advection-diffusion, are level 
curves of the 'heat function' f (Kimura & Bejan 1983), which we define in complex 
notation via Vi? = iF. For a linear diffusivity, we integrate Eq. (4.6) to obtain the 
heat function for any conformal mapping. 



H = 2'Rjef{z) 



Pe (Im/(z))erf (VP^Im/(z)) + J exp (-Pe {lmf{z)f) 



(4.7) 

which shows how the total-flux lines cross over smoothly from fluid streamlines 
outside the diffusion layer {H ^ Pe Im/(2;)^, Pe Im/(z) ^ 1) to diffusive-flux lines 
near the absorbing surface [H ~ 2-\/Pe /tt Re/(z), Pe \vtif{z) <C 1). 

On the absorbing surface, Im/(2;) = 0, the flux density is purely diffusive and 
in the normal direction. Its spatial distribution is determined geometrically by the 
conformal map, 



VeS'{Q)\f'{z)\ on lmf{z)=Q, (4.8) 

and only its magnitude depends on Pe , as predicted by the Equivalence Theorem. 
(For a linear diffusivity, 5"(0) = What appears to be the only previous result 

of this kind is due to Koplik et al. (1994, 1995) in the context of tracer dispersion 
by linear advection-diffusion in porous media. In the case of planar potential flow 
from a dipole source to an equipotential absorbing sink, they proved that the spa- 
tial distribution of surface flux is independent of Pe . Here we see that the same 
conclusion holds for all similarity solutions to Eq. (4.1), even if (i) diffusive flux 
is not directed along streamlines; (ii) the diffusivity is a nonlinear function of the 
concentration; and (iii) the domain is on a curved surface. 

(6) Streamline Coordinates 

In proving their equivalence theorem, Koplik et al. (1994, 1995) transform 
Eq. (4.1) in the linear case, 6(c) = 1, to 'streamline coordinates', 

Pe = -I- c^^, (4.9) 

where $ = + «-(/' is the complex potential, the velocity potential, and ^, the 
streamfunction. Because the independent and dependent variables are interchanged, 

t Sen & Yang (2000) have recently shown that the heat function satisfies Laplace's equation, 
V i? = 0, in certain potential-dependent coordinates, V = e~^^ V. This might seem related 
to our theorems, but it does not provide a basis for conformal mapping of the domain because the 
coordinate transformation is not analytic. Its value is also limited by the fact that the boundary 
conditions on H are not known a priori. For example, on a surface where the concentration is 
specified, the unknown flux is also required. These difficulties underscore the fact that the solutions 
of Eq. (4.1) are fundamentally non-harmonic. 
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this is a type of hodographic transformation (Whitham 1974; Ben Amar & Poire 
1999). The physical interpretation of Eq. (4.9) is that advection (the left-hand 
side) is directed along streamlines, while diffusion (the right-hand side) is also 
perpendicular to the streamlines, along iso-potential lines. In high-Reynolds-number 
fluid mechanics, this is a well known trick due to Boussinesq (1905) still in use today 
(Hunt & Eames 2002). Streamline coordinates are also used in Maksimov's method 
for dendritic solidification from a flowing melt (Cummings et al. 1999). 

Boussinesq's transformation is simply a conformal mapping to a geometry of 
uniform flow. Any obstacles in the flow are mapped to line segments (branch cuts of 
the inverse map) parallel to the streamlines. Among the solutions (4.3), streamline 
coordinates correspond to the map, f{z) = \fz, from a plane of uniform flow past 
an absorbing flat plate on the positive real axis (the branch cut), as shown in the 
top right panel of Fig. 2. In this geometry, we have the boundary- value problem, 
Pe ^ = V^c, c(x > 0, 0) = 0, c(-oo, y) = 1, which Carrier et al. (1983) have solved 
using the Wcincr-Hopf technique. More simply, Greenspan has introduced parabolic 
coordinates (as in Greenspan 1961), to immediately obtain the similarity solution 
derived above, c{x,y) = erf (v^Pe'r;), where 2rf = —x + \/ x"^ + y^. The reason 
why this solution exists, however, only becomes clear after conformal mapping to 
non-streamline coordinates in the upper half plane. (See also Cummings et al. 1999.) 

As this example illustrates, streamline coordinates are not always convenient, 
so it is \iscful to exploit the possibility of conformal mapping to other geometries. 
For similarity solutions, it is easier to work in a plane where the diffusive flux 
lines are parallel. Streamline coordinates are also often poorly suited for numerical 
methods because stagnation points are associated with branch-point singularities. 
This is especially problematic for free boundary problems: For flows toward infinite 
dendrites, it is easier to determine the evolving map from a half plane (Cummings et 
al. 1999); for flows past finite growing objects, it is easier to map from the exterior 
of the unit circle (Bazant, Choi & Davidovitch 2003). 

(c) Non-similarity Solutions for Finite Absorbing Objects 

It is tempting to try to eliminate the plate from the cylindrical rim in Fig. 3 

by conformal mapping from the exterior of a finite object to the upper half plane. 
Any such mapping in Eq. (4.3), however, requires a quadrupole point source of 
flow (mapped to oo) on the object's surface. This is illustrated in the lower right 
panel of Fig. 2 by a Mobius transformation from the exterior of the unit circle, 
f{z) = (1 -|- z)/i{l — z), where a source at z = 1 ejects concentrated fluid in the 
+1 direction and sucks in fluid along the ±i directions. Thus we see that, due to 
the boundary conditions at oo, uniform flow past an absorbing cylinder (or any 
other finite object) is in a different class of solutions, where the diffusive flux lines 
depend nontrivially on Pe . In streamline coordinates, this includes the problem of 
uniform flow past a finite absorbing strip, which requires solving Wijngaarden's 
integral equation (Cummings et al. 1999). 

Here, we study only the high-Pe asymptotics of advection-diffusion layers around 
finite absorbing objects. Consider again the example of fiow past a cylindrical rim on 
a flat plate (Fig. 3). Because disturbances in the concentration decay exponentially 
upstream beyond an 0(Pe ~^/^) distance, removing the plate on the downstream 
side of the cylinder has no effect in the limit Pe ^ oo, except on the plate itself 
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(the branch cut), so the solution (4.4)-(4.5) is also asymptotically valid near a 
finite absorbing cylinder (without the plate). 

More generally, if z = h{q) is the conformal map from the exterior of any singly 
connected finite object to the exterior of the unit circle, then the non-harmonic 
concentration field has the asymptotic form, 



c(g, q) ~ erf 



Pe Im 7%) 



(4.10) 



as Pe — > oo everywhere except in the wake near the pre-image of the positive real 
axis, a branch cut corresponding to the 'false plate'. The convergence is not uniform, 
since the false plate always spoils the approximation sufficiently far downstream, for 
a fixed Pe ^ 1. The validity of Eq. (4.10) near the surface of the object, however, 
allows us to calculate the normal flux density using Eq. (4.8), 




n-Vc~2J — sin - (4.11) 



as Pe oo for all 9 = argh{q) ^ Pc^^^^ away from the rear stagnation point, 
= 0. The limiting Nussclt number, Nu ~ 8-\/Pe /tt, is also easily calculated by 
mapping the rim (with the false plate) to the upper half plane where the normal 
flux density is uniform, 2-y/Pc /tt, on a line segment of length four (from -2 to 2). 

As explained in section 2(c), Equation (4.11) describes the non-harmonic prob- 
ability measure for fractal growth by steady advection-diffusion in a uniform po- 
tential flow in the limit Pe ^ oo. This model, which we might call 'advection- 
diffusion- limited aggregation' (ADLA), is perhaps the simplest generalization of 
the famous DLA model of Witten and Sander (1981) allowing for more than one 
bulk transport process. The resulting competition between advection and diffusion 
produces a crossover between two distinct statistical 'phases' of growth. As expected 
from renormalization-group theory (Goldenfeld 1992), the crossover connects 'fixed 
points' of the growth measure, describing self-similar dynamics. For small initial 
Peclet numbers, Pe (0) ^ 1, the growth measure of ADLA is well approximated 
by the uniform harmonic measure of DLA and the concentration by the similar- 
ity solution, c{q,q) oc lm\ogh{q), but this is an unstable fixed point. Regardless 
of the initial conditions, the Peclet number diverges, Pe (t) = UL{t)/D oo, as 
the object grows, so the concentration eventually approaches the new similarity 
solution in Eq. (4.10). At this advection-dominated stable fixed point, the growth 
measure obeys Eq. (4.11). The sin 0/2 dependence causes anisotropic fractal growth 
at long times favoring the direction of incoming, concentrated fiuid, 9 = w, and the 
total growth rate (Nu ) is proportional to ^Pc (t). Such analytical results serve to 
illustrate the power of conformal mapping applied to systems of invariant equations. 



5. Electrochemical Transport 

(a) Simple Approximations and Conformal Mapping 

Conservation laws for gradient-driven fluxes also describe ionic transport in di- 
lute electrolytes. Because the complete set of equations and boundary conditions 
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(below) are nonlinear and rather complicated, the classical theory of electrochem- 
ical systems involves a hierarchy of approximations (Newman 1991). Conformal 
mapping has long been applied in the simplest case where the current density, J, is 
proportional to the gradient of a harmonic function, </>, the electrostatic potential 
(Moulton 1905; Hine 1956). 

This approximation, the 'primary current distribution', describes the linear re- 
sponse of a homogeneous electrolyte to a small applied voltage or current, as well 
as more general conduction in a supporting electrolyte (a great excess of inactive 
ions). The assumptions of Ohm's Law, J = crE = — aVi/) (with a constant conduc- 
tivity, a) and no bulk charge sources or sinks, V • J = 0, are analogous to those of 
potential flow and incompressibility describe above. Each electrode is assumed to 
be an equipotential surface (sec below), so the potential is simply that of a capaci- 
tor — harmonic with Dirichlet boundary conditions. Naturally, classical conformal 
mapping from electrostatics (Churchill & Brown 1990; Needham 1997) have been 
routinely applied, but it seems conformal mapping has never been applied to any 
more realistic models of electrochemical systems. 

The 'secondary current distribution' introduces a kinetic boundary condition, 
fb -J = R{4>), which equates the normal current with a potential-dependent reaction 
rate, e.g. given by the Butler- Volmer equation (see below). In this case, conformal 
mapping could be of some use. Although the boundary condition acquires a non- 
constant coefficient, |/'|, from Eq. (2.7), Laplace's equation is preserved. 

A more serious complication in the 'tertiary current distribution' is to allow 
the bulk ionic concentrations to vary in space (but not time). Ohm's law is then 
replaced by a nonlinear current-voltage relation. Our main insight here is that 
conformal mapping can still be applied in the usual way, even though the equations 
are nonlinear and the potential, non-harmonic. 

(&) Dilute-Solution Theory 

In the usual case of a dilute electrolyte, the ionic concentrations, {ci, C2, . . . , Cjv}, 
and the electrostatic potential, (p, satisfy the Nernst-Planck equations (Newman 
1991), which have the form of Eqs. (3.1) and (3.2), where the 'advection' velocities, 
Uj = —ZiCfii^cj), are due to migration in the electric field, E = — V(/). Here, Zje 
is the charge (positive or negative) and /Uj the mobility of the ith ionic species. 
The diffusivitics are given by the Einstein relation, _D, = fcsT/ij, where k-B is 
Boltzmann's constant and T, the temperature. Scaling concentrations to a reference 
value, C, potential to the thermal voltage, kT/e, length to a typical electrode 
separation, L, and assuming that Di, T, and e are constants, the steady-state 
equations take the dimensionlcss form, 

V^ci + V • (ci V0) = 0. (5.1) 

The ionic flux densities are = — Vc^ — z^qV^ (scaled to DiC/L). 

Because dissolved ions are very effective at charge screening, significant diffuse 
charge can only exist in very thin (1 — lOOnm) interfacial double layers, where 
boundary conditions break the symmetry between opposite charge carriers. The 
'bulk' potential (outside the double layers) is then determined implicitly by the 
condition of electroneutrality (Newman 1991), X^^Li ^^iCj = 0, which is trivially 
conformally invariant. Therefore, the most common model of steady electrochemical 
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transport, Eq. (5.1), satisfies the assumptions of the Conformal Mapping Theorem 
for any number of ionic species {N >2). Although the equations differ from those of 
advection-diffusion in a potential flow, wc can still map to electric- field coordinates 
(the analog of streamline coordinates), or any other convenient geometry. 

Although the equations are conformally invariant, the boundary conditions are 
so only in certain limits. General boundary conditions express mass conservation, 
either n • Fj = for an inert species, or 

n- VFi =i?i(ci,(/)), (5.2) 

for an active species at an electrode, where Ri{ci,(j)) is the Faradaic reaction-rate 
density (scaled to DiC/L). It is common to assume Ahrrcnius kinetics, 

Riia, (p) = fc+Cie^'«+("^-'^=) - k-CrS-^*"-^*-*-^ (5.3) 

where fc+ and /c_ arc rate constants for deposition and dissolution, respectively 
(scaled to Di/L), a± are transfer coefficients, is the concentration of the reduced 
species (scaled to C) and (pe is the electrode potential (scaled to kT/e). Taking 
diffuse interfacial charge into account somewhat modifies R(ci,(f>), but the basic 
structure of Eq. (5.2) is unchanged (Newman 1991; Bonnefont et al. 2001). Confor- 
mal mapping introduces a non-constant coefficient, |/'|, in Eq. (5.2), but conformal 
invariance is restored in the case of 'fast reactions' (A:_|_ ^ l,fc_Cr 1), in which 
equilibrium conditions prevail, i? = 0, even during the passage of current. For a 
single active species (say i = l), the bulk potential at an electrode is then given by 
the (dimensionless) Nernst equation. 



^4>e, = - ^,"'\ : „ (5.4) 



where k = k+/k-Cr is an equilibrium constantf. 



(c) Conformal Mapping with Concentration Polarization 

The voltage across an electrochemical cell is conceptually divided into three 
parts (Newman 1991): (i) the 'Ohmic polarization' of the primary current distribu- 
tion, (ii) the 'surface polarization' of the secondary current distribution, and (iii) 
'concentration polarization', the remaining voltage attributed to non- uniform bulk 
concentrations. Although concentration polarization can be significant, especially 
at large currents in binary electrolytes, it is difficult to calculate. Analytical results 
are available only for very simple geometries (mainly in one dimension), so our 
method easily produces new results. 

For example, consider a symmetric binary electrolyte {N = 2) of charge number, 
z = Z-^- = —Z-, where the concentration, c — c+ = c_, and the potential satisfy, 

V^c = and V • (cV^i) = 0. (5.5) 

(The concentrations are harmonic only for A'^ = 2.) Assuming that anions are 
chemically inert yields an invariant zero-flux condition at each electrode, h- (cV(/> — 

t Expressing Eq. (5.3) in terms of the 'surface overpotential' , rjs = (p — 4'e — yields the 

more famihar Butler- Volmer equation (Newman 1991). 
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Figure 4. Similarity solutions for the electrostatic potential (contour plot) and cur- 
rent/electric-field lines (yellow) in a binary electrolyte at 90% of the limiting current 
{k — 1). The simple solution for parallel-plate electrodes (top) is conformally mapped to 
semi-infinite plates (bottom left) and misaligned coaxial cylinders (bottom right). 

Vc) — 0, and (to break degeneracy) a constraint on the integral of c, which sets the 
total number of anions (Bonnefont et al. 2001). In the limit of fast reactions, the 
bulk potential at each electrode is given by the Nernst equation, cj) — (p^ ~ log fcc, 
where we scale 4> to ksT / ze and assume a+ — a_ = 1. 

A class of similarity solutions is obtained by conformal mapping, w = f{z), to 
a strip, —1 < Imui < 1, representing parallel-plate electrodes. We set 0e = at the 
cathode (Imw — —1) and 0e = V, the applied voltage (in units of ksT/ze), at the 
anode (Imw — 1). We then solve c" = and (c</>')' = with appropriate boundary 
conditions to obtain a general solution for any conformal mapping to the strip: 

c=l + JIm/(z), 0-log( ^ + ^:^y ), (5.6) 

where J = tanhiV / 4) is the uniform current density in the strip, scaled to its 
limiting value, Jum — 2zeD+C / L. As J ^ 1, strong concentration polarization 
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develops near the cathode, as shown in Fig. 4 for J = 0.9. At J = 1, the bulk 
concentration at the cathode vanishes, and the cell voltage diverges due to diffusion 
limitation. 

The classical conformal map, z = f~^{w) = nw+e^^ (Churchill & Brown 1990), 
unfolds the strip like a 'fan' to cover the 2;-plane and maps the electrodes onto two 
half plates (Imz = ±7r , Rcz < —1). As shown in Fig. 4, this solution describes the 
fringe fields of semi-infinite, parallel-plate electrodes. The field and current lines are 
cycloids, Za{ri) = 7ra + iwr] + e^°'e™'^, as in the limit of a harmonic potential at low 
currents, (j) ^ J Im f{z) — \ogk. At high currents, the magnitude of the electric field 
is greatly amplified near the cathode (the lower plate) by concentration polarization, 
but the shape of the field lines is always the same. This conclusion also holds for 
all other conformal mappings to the strip, such as the Mobius-log transformation, 
w = f{z) = i(l -|- log(5z — 3)/(5 — 3z)), in Fig. 4 from the region between two 
non-concentric circles. 

It is interesting to note that the Equivalence Theorem applies to some physi- 
cal situations and not others. Similarity solutions like the ones above can only be 
derived for two equipotential electrodes by conformal mapping to a strip, where 
the current is uniform. In all such geometries, the electric field lines have the same 
shape as in the primary current distribution. For three or more equipotential elec- 
trodes, however, this is no longer true because conformal mapping to the strip is 
topologically impossible, and thus similarity solutions do not exist. When the bulk 
potential varies at the electrodes according to Eq. (5.2), the electric field lines gen- 
erally differ from both the primary and secondary current distributions, even for 
just two electrodes. 

6. Conclusion 

We have observed that the nonlinear system of equations (2.4) involving 'dot prod- 
ucts of two gradients' is conformally invariant. This has allowed us to extend the 
classical technique; of conformal mapping to some non-harmonic functions arising in 
physics. Examples from transport theory are steady conservation laws for gradient- 
driven fluxes, Eq. (3.2). For one variable, the equations in our class (including some 
familiar examples in nonlinear diffusion) can always be reduced to Laplace's equa- 
tion. For two or more variables, the general solutions are not simply related to 
harmonic functions, but all similarity solutions exhibit an interesting geometrical 
equivalence. 

For two variables, there is one example in our class, steady advection-diffusion 
in a potential flow, to which conformal mapping has previously been applied. In this 
case, our method is equivalent to Boussinesq's streamline coordinates, but some- 
what more general. A nonlinear diffusivity is also allowed, and the mapping need 
not be to a plane of uniform flow (parallel streamlines). In a series of examples, 
we have considered flows past absorbing leading edges and have generalized a re- 
cent equivalence theorem of Koplik, Redner, and Hinch (1994, 1995). We have also 
considered the flows past finite absorbing objects at high Peclet number. 

Our class also contains the Nernst-Planck equations for steady, bulk electro- 
chemical transport, for which very few exact solutions are known in more than one 
dimension. In electrochemistry, conformal mapping has been applied only to har- 
monic functions, so we have presented some new results, such as the concentration 
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polarizations for semi- infinite, parallel-plate electrodes and for misaligned coaxial 
electrodes. More generally, we have shown that Ohm's Law gives the correct spatial 
distribution (but not the correct magnitude) of the electric field on any pair of 
equipotential electrodes in two dimensions, even if the transport is nonlinear and 
non-Laplacian, although this is not true for three or more electrodes. Such results 
could be useful in modeling micro-electrochemical systems, where steady states arc 
easily attained (due to short diffusion lengths) and quasi-planar geometries are 
often arise. 

As mentioned thoughout the paper, our results can be applied to a broad class 
of moving free boundary problems for systems of non-Laplacian transport equations 
(Bazant, Choi & Davidovitch 2003). In contrast, the vast literature on conformal- 
map dynamics (cited in section 1) relies on complex-potential theory, which only ap- 
plies to Laplacian transport processes. Nevertheless, standard formulations, such as 
the Polubarinova-Galin equation for continuous Laplacian growth (Howison 1992) 
and the Hastings-Lcvitov (1998) method of iterated maps for DLA, can be eas- 
ily generalized for coupled non-Laplacian transport processes in our class. In the 
stochastic case, non-harmonic probability measures for fractal growth can be de- 
fined on any convenient contour, such as the unit circle. As an example, wc have 
derived the stable fixed point of the growth measure for an arbitrary absorbing 
object in a uniform background potential flow, Eq. (4.11). This sets the stage for 
conformal-mapping simulations of ADLA, which might otherwise seem intractable. 
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